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Abstract 

We consider numerically the depinning transition in the random- 
field Ising model. Our analysis reveals that the three and four di- 
mensional model displays a simple scaling behavior whereas the five 
dimensional scaling behavior is affected by logarithmic corrections. 
This suggests that d = 5 is the upper critical dimension of the de- 
pinning transition in the random-field Ising model. Furthermore, we 
investigate the so-called creep regime (small driving fields and tem- 
peratures) where the interface velocity is given by an Arrhenius law. 

Keywords: Disordered media, Driven interfaces, Numerical investigations, 
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1. INTRODUCTION 

Using the random-field Ising model (RFIM) we have studied the dynamics 
of driven interfaces in disordered ferromagnetic systems by means of Monte 
Carlo simulations. The interface separates regions of opposite magnetization 
and it is driven by a homogeneous field H. Without thermal fluctuations 
(T = 0) its dynamics is affected by the disorder in the way that the interface 
moves only if a sufficiently large driving field is applied. With decreasing 
driving field the interface velocity vanishes at a critical threshold H c and 
is pinned below H c , due to the random- field. The field dependence of the 
interface velocity is given by 

v(h) ~ h p (1) 

with h = (H — H c )/H c > 0, i.e. this pinning/ depinning transition can be 
considered as a continuous phase transition. A sketch of this behavior is 
shown in Fig. [I]. 
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The paradigm of the depinning transition is the Edwards- Wilkinson equa- 
tion ( [Edwards and Wilkinson] , |1982| ) with quenched disorder (QEW). The 
critical exponent (3 of the QEW equation has been estimated within an renor- 
malization group approach (|Nattermann, Stepanow, Tang and Leschhorn|, 



|1992j ) It is expected that the critical behavior of both the RFIM and the 
QEW equation is characterized by the same exponents, i.e. both models are 
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Barabasi and Stanley! , |1994| ). 

In this paper we determine numerically the critical exponent (3 of the 
depinning transition in the RFIM in various dimensions. We compare our 
results with those of the renormalization group approach mentioned above. 
In particular, we consider the five dimensional case where logarithmic correc- 
tions to the above scaling ansatz, Eq. ([[]), occur. These logarithmic correc- 
tions are a typical feature of the upper critical dimension where the so-called 
mean field solution becomes valid. 

For finite temperatures no pinning of the interface occurs due to the 
thermal fluctuations. Even for driving fields well below H c the interface is 
expected to exhibit a so-called creep motion, see Fig. [I], in which the interface 
velocity obeys an Arrhenius law. This Arrhenius law has been found for the 
QEW equation flChauve, Giamarchi and Doussal , 1998|) and we compare 
these results with those of the RFIM. 



2. THE RFIM MODEL 

In order to study the dynamics of moving interfaces numerically, we have 
investigated the rf-dimensional RFIM which is characterized by the Hamilto- 
nian, 

n = - J - Y / s t .s j -HY / s l -Y,h l s l . (2) 

(id) 

The sum in the first term is taken over all pairs of neighbored spins. It de- 
scribes the exchange interaction between neighboring spins (Si = ±1) and a 
parallel alignment of nearest neighbors is energetically favored. Additionally, 
the spins are coupled to a homogeneous driving field H and to a quenched 
local random-field hi with (hihj) oc 5y and (hi) = 0. The probability density 
p that the local random-field takes some value hi is given by 

v(h)-{ (2A)_1 f ° r N<A (3) 
P^)-\ o otherwise, [6) 

i.e. hi is uniformly distributed. Throughout the paper we consider the case 
A = 1.7. Using both simple cubic (sc) and body centered cubic (bcc) lat- 
tices with antiperiodic boundary conditions we started each simulation for 
an initially flat interface. In non-disordered systems (A = 0) the interface 
moves for any finite driving field. This limiting behavior for A — > can be 
recovered if the interface moves along the diagonal direction of simple cubic 
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(sc) lattices or along the i-direction of body centered cubic (bcc) lattices, for 
details see flNowak and Usadel| , |199S| ). 

In our Monte Carlo simulations the dynamics of the interface motion is 
given by a Glauber dynamics with transition probabilities according to a 
heat-bath algorithm, see e.g. (Pinder and Heermann| , |1997|) . The interface 
velocity is the basic quantity in our investigations. Since a moving interface 
corresponds to a magnetization M which increases with time t (given in 
Monte Carlo steps per spin) the interface velocity is obtained according to 
v ~ (dM/dt), where (...) denotes an appropriate disorder average. 



3. CRITICAL INTERFACE DYNAMICS 

In the following we analyse the interface dynamics in the vicinity of the 
critical threshold H c at zero temperature (T = 0) on a bcc lattice. The 
interface velocity as a function of the driving field for the three dimensional 
model is shown in Fig. 0. The data suggest that v vanishes continuously 
at the critical threshold H c . Fitting the data to the above ansatz, Eq. ([[]), 
one obtains (3 = 0.642 ± 0.06 and H c = 1.357 ± 0.008. This agrees with 
previous results. In a simple cubic lattice an exponent (3 = 0.6 ± 0.11 was 
observed for an interface moving along the [001] direction ( [Amaral, Barabasl 



|and Stanley! , |1994]) , and (3 = 0.66 ±0.04 was observed for an interface moving 
along the [111] direction ( [Roters et gl.[ |1999|) , respectively. Additionally, 
these values agree within the error-bars with the value (3 = 2/3 found for the 
QEW equation with quenched disorder ( [Nattermann, Stepanow, Tang and 
ILeschhomj gggg ). 



The above scaling behavior of the interface velocity, Eq. (JJ), is also found 
in two ( pNowak and Usadel| , |1998| ) and four dimensions (not shown). But 
in five dimensions the simple scaling ansatz Eq. ([[]) does not describe the 
velocity-field dependence in a convincing manner. In the inset of Fig. |3] we 
show the logarithmic derivative /3 e g = dlnv/dlnh for various driving fields. 
Instead of a saturation for h — ► the effective exponent (3 e g displays a strong 
curvature which contradicts the scaling ansatz in Eq. ([!]). 

In the following we will show that the data for d = 5 are consistent with 
the assumption that the scaling behavior is given by Eq. ([!]) modified by 
logarithmic corrections. Usually logarithmic corrections affect the scaling 
behavior at the upper critical dimension d c . Below d c the critical exponents 
depend on the dimension of the system while for d > d c they are independent 
of the particular dimension and are given by those of the mean field theory. 
At d c itself the mean field scaling behavior is modified by logarithmic cor- 
rections. The above mentioned renormalization group approach to the QEW 
equation predicts the values d c = 5 and Pmf = 1- 

Inspired by the usual scaling behavior at the critical dimension obtained 
from renormalization group approaches, see e.g. ( [Brezin, Le Guillou and 



Zinn- Justing p.976|) , we assume that the asymptotic (h — > 0) scaling behavior 
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of the velocity is given in leading order by 

v ~ \h\^F\\nh\ Xh . (4) 

with /?mf — 1 and some unknown exponent Xh- The corresponding data as 
well as the fit according to Eq. (^) are shown in Fig. [3| From the fit which 
works for small values of h we obtain H c = 1.143 ± 0.01 and Xh = 0.42 ±0.1. 
It remains to be shown that the scaling behavior above d c is given by the 
pure mean field scaling behavior. Simulations of the six dimensional model 
are in progress and the results will be published elsewhere. 



4. THE CREEP MOTION 

In the previous section we focused our attention to the critical behavior at 
zero temperature. Including thermal fluctuations no pinning of the interface 
takes place, see Fig. |T[ For example, at H = H c the velocity grows with 
the temperature according to v ~ T l ' s ( |Nowak and Usadelj , [1998| ; pHir" 



et al.[ |1999| ). The interface moves even for driving fields well below the 
critical threshold. In this case the interface velocity is expected to obey an 
Arrhenius law 

V - H //•'/') exp I t(H) 



T 



(5) 



which is characterized by its prefactor C(H,T) and an effective energy bar- 
rier E(H). Recently flChauve, Giamarchi and Doussa| , |1998|) investigated 



the creep regime in the QEW equation and found that the prefactor of the 
Arrhenius law is given by 

C{H,T) =c{H)T~ x (6) 

with c(H) being some function of the driving field. The temperature de- 
pendence is characterized by an exponent x. We apply this ansatz to our 
numerically obtained interface velocities in the 3d RFIM. Figure |] shows the 
data rescaled according to Eq. @ and (||). As can be seen from the data, we 
obtain nearly straight lines for x = 0.79 ± 0.09 suggesting that the tempera- 
ture dependence of the prefactor in Eq. (|5D as well as the Arrhenius law itself 
are a proper description of the velocity in the creep regime. The slope of 
the lines in Fig. f| corresponds to the effective energy barrier E(H). Plotting 
\nvt x vs. E(H)/T the data collapse onto a single curve, see Fig. ||, without 
making any assumption about c(H). This suggests that c(H) »s const holds. 
However, details of our results do not coincide with those of QChauve 



|Giamarchi and Doussal , 1998| ). In particular, these authors observed a sig- 



nificant field dependence of the prefactor c(H) ~ with /i > and found 
a negative exponent x. Especially the opposite sign of x is remarkable. 
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5. SUMMARY 

We investigated the dynamics of driven interfaces for the random-field Ising 
model (RFIM). In our analyses we focus on the field and temperature de- 
pendence of the interface velocity and we consider the depinning transition 
especially at temperature T = 0. In the three dimensional model we find 
that v(H) obeys a power law characterized by an exponent (5. In d = 5 the 
scaling behavior can be described by the mean field exponent modified by 
logarithmic corrections. This suggests that the upper critical dimension of 
the depinning transition in the RFIM is d c = 5. 

We also investigate the interface dynamics in the so-called creep regime 
which occurs for driving fields well below the critical threshold and small but 
finite temperatures. As predicted by a renormalization group approach the 
creep behavior is characterized by an Arrhenius law. But the details of the 
Arrhenius law differ from the predicted behavior. 
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Figure 1: Sketch of the interface velocity v on the driving field H and 
temperature T. The bold line corresponds to T = 0. With increasing tem- 
perature the depinning transition is smeared out (thin solid line). In the 
creep regime the interface velocity obeys an Arrhenius law. 




Figure 2: Dependence of the interface velocity v on the driving field H in 
the three dimensional model for different system sizes L. The data can be 
fitted (dashed line) according to Eq. ([I]). From the fit one obtains H c = 
1.357 ± 0.008 and (3 = 0.642 ± 0.06. 
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Figure 3: Dependence of the interface velocity v on the driving field H in 
the five dimensional model for different system sizes L. The dashed line is 
a fit according to Eq. where /?mf = 1 has been used. As can be seen 
from the data, the fit works quite well for H —>■ H c and one obtains H c = 
1.1425±0.001 and Xh = 0.42±0.1. The inset shows the logarithmic derivative 
/3 e ff = d\nv/d\nh. The strong curvature indicates that the velocity field 
dependence cannot be described by Eq. ([[]). 




Figure 4: Interface velocities obtained at different driving fields H = 
0.3, 0.35, 0.4, 0.6 and rescaled according to Eq. (J|) and (||). Varying x 
we obtain nearly straight lines for x = 0.79 ± 0.09. 
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Figure 5: Rescaled interface velocities. The data coincide with those shown 
in Fig. |j. They are rescaled with the numerically determined energy barrier 
E(H). The coalescence of the data suggests c(H) « const. 



